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Keno Eilers (-^^E Betti Hartmann (''^Q Valeria Kagramanova ("'0 Isabell Schaffer C*)! and Catalin Toma (^Q 

(a) Institut fiir Physik, Universitdt Oldenburg, 26111 Oldenburg, Germany 

(b) School of Engineering and Science, Jacobs University Bremen, 28759 Bremen, Germany 

(Dated: April 23, 2013) 

We study the geodesic motion of test particles in the space-time of non-compact boson stars. 
These objects are made of a self-interacting scalar field and - depending on the scalar field's mass 
- can be as dense as neutron stars or even black holes. In contrast to the former these objects do 
not contain a well-defined surface, while in contrast to the latter the space-time of boson stars is 
globally regular, can - however - only be given numerically. Hence, the geodesic equation also has 
to be studied numerically. We discuss the possible orbits for massive and massless test particles 
and classify them according to the particle's energy and angular momentum. The space-time of a 
boson star approaches the Schwarzschild space-time asymptotically, however deviates strongly from 
it close to the center of the star. As a consequence, we find additional bound orbits of massive test 
particles close to the center of the star that are not present in the Schwarzschild case. Our results 
can be used to make predictions about extreme-mass-ratio inspirals (EMRIs) and we hence compare 
our results to recent observational data of the stars orbiting Sagittarius A* - the radiosource at the 
center of our own galaxy. 
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I. INTRODUCTION 

Solitons play an important role in many areas of physics. As classical solutions of non-linear field theories, they are 
localized structures with finite energy, which are globally regular. In general, one can distinguish between topological 
and non-topological solitons. While topological solitons [l| possess a conserved quantity, the topological charge, that 
stems (in most cases) from the spontaneous symmetry breaking of the theory, non-topological solitons [2|, |3| have a 
conserved Noether charge that results from a symmetry of the Lagrangian. The standard example of non-topological 
solitons are Q-balls [^], which are solutions of theories with self- interacting complex scalar fields. These objects are 
stationary with an explicitly time-dependent phase. The conserved Noether charge Q is then related to the global 
phase invariance of the theory and is directly proportional to the frequency. Q can e.g. be interpreted as particle 
number [2]. While in standard scalar field theories, it was shown that a non-renormalizable ^^-potential is necessary 
[31 , supersymmetric extensions of the Standard Model (SM) also possess Q-ball solutions [0| . In the latter case, several 
scalar fields interact via complicated potentials. It was shown that cubic interaction terms that result from Yukawa 
couplings in the superpotential and supersymmetry (SUSY) breaking terms lead to the existence of Q-balls with non- 
vanishing baryon or lepton number or electric charge. These supersymmetric Q-balls have been considered aspossible 
candidates for baryonic dark matter [3] and their astrophysical implications have been discussed Q- In 0, these 
objects have been constructed numerically using the exact form of a scalar potential that results from gauge-mediated 
SUSY breaking. However, this potential is non-differcntiable at the SUSY breaking scale. In [l^ a differentiable 
approximation of this potential was suggested and the properties of the corresponding Q-balls have been investigated 
in 3 + 1 dimensions. This was extended to d -I- 1 dimensions in [ll|. Q-ball solutions with a $^-potential in 3 -I- 1 
dimensions have been studied in detail in [3, [l^, Il3 ■ 

It was realized that next to non-spinning Q-balls, which are spherically symmetric, spinning solutions exist. These 
are axially symmetric with energy density of toroidal shape and angular momentum J = ttiQ, where Q is the Noether 
charge of the solution and to G Z corresponds to the winding around the z-axis. Approximated solutions of the non- 
linear partial differential equations were constructed in |5j| by means of a truncated series in the spherical harmonics 
to describe the angular part of the solutions. The full partial differential equation was solved numerically in |12h14| . 
It was also realized in Q that in each m-sector, parity-even (P = -1-1) and parity-odd {P = —1) solutions exist. 
Parity-even and parity-odd refers to the fact that the solution is symmetric and anti-symmetric, respectively with 
respect to a refiection through the x-y-plane, i.e. under ^ it — 9. Complex scalar field models coupled to gravity 
possess so-called "boson star" solutions |15l - [2ll | . In [l^, [Ij, [24I boson stars have been considered that have flat space- 
time limits in the form of Q-balls. These boson stars are hence self-gravitating Q-balls. In [ll|, [2j] the gravitating 
generalizations of the supersymmetric Q-balls studied in [10| have been discussed in d -I- 1 dimensions. It was found 
that the behaviour of the mass and charge at the critical value of the frequency depends crucially on the number of 
dimensions d. While in most models considered, the scalar field function is exponentially decaying and hence different 
notions of a boson star radius exist, this is different in models with a V-shaped potential [2J,[2a]- In this case, compact 
boson stars with a well-defined outer radius (very similar to those of "standard stars" ) can be given [26 28] . 

In this paper we are interested in the possibility to detect boson stars through the motion of massive and massless 
test particles in their space-time. In particular, we will be interested in the difference between test particle motion 
in a boson star space-time and a Schwarzschild space-time. Since boson stars are very compact, they have been 
considered as alternatives to supermassive black holes [29| residing e.g. in the center of galaxies and geodesic motion 
of massive test particles describing extreme-mass-ratio inspirals (EMRIs) has been discussed in some particular cases 
in [3^. While objects with a well-defined surface as alternatives to the supermassive black hole at the center of our 
own galaxy, the Milky Way have been ruled out [31] , the boson stars studied in this paper have a scalar field falling of 
exponentially at infinity and hence strictly speaking do not have a well-defined outer surface outside which the energy 
density and pressure, respectively, vanishes. We hence study non-compact boson stars in this paper and we will make 
a detailed analysis of the motion of massive and massless test particles in the space-time of such a boson star. 

Our paper is organised as follows: in Section II, we give the field theoretical model, the Ansatz and the equations of 
motion to describe the space-time of a non-spinning, non-compact boson star. In Section III we discuss the geodesic 
equation and give our numerical results in Section IV. We conclude in Section V. 

II. THE SPACE-TIME OF A NON-COMPACT BOSON STAR 

In the following we will discuss the field theoretical model to describe the space-time of a non-spinning, non-compact 
boson star in which the test particles will move on geodesies. 
The action S of the field theoretical model reads: 



where R is the Ricci scalar, G denotes Newton's constant and the matter Lagrangian is given by 

£,„ = -d^^d^^* - Um) (2) 

where $ denotes a complex scalar field and we choose as signature of the metric ( — h ++)• C^d'l'l) is the scalar field 
potential that arises in gauge-mediated supersymmetric breaking in the Minimal Supersymmetric extension of the 
Standard Model (MSSM) ^^: 

um) ^mWsusy(^~^^p(-^)) , (3) 

\ \ '/susy / / 

where rysusy corresponds to the scale below which supersymmetry is broken and is roughly 1 TeV, while m denotes 
the scalar boson mass. For |<l>p ^ 'ysusy '^^ '^^^ Si^e a polynomial expansion of this potential which reads 

um) = mm' - ^m' + ^m' + om') ■ w 

/susy /susy 

Neglecting terms of |$|^ and higher this potential has the form of the scalar potential frequently used in the study of 
Q-balls and boson stars [J, la| and reads 

?7(|$|) = A|$p(7|0|^-a|$|2+/3) , (5) 



where A, 7, a and /3 are constants such that \/A/3 corresponds to the scalar boson star mass m. In previous studies, 
the constants have been fixed to particular values which fulfill the requirements to find non-trivial solutions to the 
equations of motion. 

The matter Lagrangian Cm © is invariant under the global U(l) transformation 

$ -^ ^e'^^ , X constant . (6) 

As such the locally conserved Noether current j^, /i = 0, 1, 2, 3, associated to this symmetry is given by 

f = -i {^*^^''^ - $a^$*) with j^l = . (7) 

The globally conserved Noether charge Q of the system then reads 

Q^~ j ^—gfd\. (8) 

Finally, the energy-momentum tensor is given by 



Tfj.„ — g^i,C 2^ 



-Sm-^ 



]^g"p {dmdp^ -f dp^*dm + um 



dmdu<^ + dmdp^ . (9) 



The coupled system of ordinary differential equations is given by the Einstein equations 

G^, = 8^GT^, (10) 

with T^i, given by ([9]) and the Klein-Gordon equation 

In the following, we want to study non-spinning boson stars. For the metric we use the following Ansatz in isotropic 
coordinates 

ds^^~f{r)de+^-P-Adr'+r'de^+r'^i^'Sd^'] ■ (12) 

/('') 

For the complex scalar field, we use a stationary Ansatz that contains a periodic dependence of the time-coordinate t: 

a>(i,r) = e^"*0(r), (13) 



where a; is a constant and denotes the frequency. 

In order to be able to use dimensionless quantities we introduce the following rescalings 



r 
m 



muj , (f)^ risnsy4> 



and find that the equations depend only on the dimensionless coupling constant 

„2 



K 



SirGrii. 



Stt. 



Mp^i 



(14) 



(15) 



where Mpi is the Planck mass. Note that k denotes the ratio between the SUSY breaking scale and the Planck mass 
which we expect to be small k ^ lO^'^^ if the SUSY breaking scale is on the order of 1 TeV. In this paper, we will 
however also study larger values of k in order to understand the qualitative behaviour of the solutions in curved 
space-time. Note that with these rescalings the scalar boson mass m becomes equal to unity. 
Using the rescalings the Einstein equations read 



2fl 



AKff (1 - exp(-^2)^ _ ^nuj'^fct? + -flf - 21 f 



fl'f 



and 



1 1 
2/1 



8kI^ (1 - exp(-02)) _ 8kw2^02 _^ ^fiif _ j;/2 



/ 



while the Klein-Gordon equation for the scalar field function is 



1 1 
2/1 



2^2 (/) - 2Z2(/)exp(-02) -t- -m' + //>' 
/ -r 



(16) 



(17) 



(18) 



The prime here and in the following denotes the derivative with respect to r. 

The equations (IT5|) - (IT51) have to be solved subject to appropriate boundary conditions. These are given by the 
requirement of regularity at the origin 



/'|r=0 = , l'\r=^=^ , </)'! 

and by the requirement of finite energy, asymptotically fiat solutions 

/(r = oo) = 1 , l{r = cx)) ~ 1 , (/) 



.0-0 



oo 







(19) 



(20) 



The set of coupled, nonlinear ordinary differential equations (|16l) - p8)) has been studied in detail in [23|, however in 
Schwarzschild-like coordinates. We will quote the relevant results in the following when necessary, but refer the reader 
to that paper for details on the field theoretical solutions. 

The behaviour of the metric function / determines the mass M of the solution. This is given by [IJ, ll^l 



M = i hm r^drf , 

Z r—^oo 



while the exphcit expression for the Noether charge reads 



Vl 2 ,2 



/ 



Note that the mass M is measured in units of AA/m. 

At spatial infinity r -^ oo the scalar field function decays exponentially 



(21) 



(22) 



(r >> 1) = - exp(-Vl -cj^r) + 



(23) 



Hence (/)(r) ^ for r < c». As such the boson star does not possess a well-defined surface and different notions of 
the radius of the boson star are possible (see discussion below). While boson stars with well-defined surface do exist 
[28| (so-called "compact boson stars"), which hence rather resemble astrophysical objects, boson stars without proper 
surface could play an important role as alternatives to supermassive black holes in the center of galaxies. 



III. THE GEODESIC EQUATION 

We consider the geodesic equation 



dxd.^^^ (24) 



dT2 P" dT dT 

where Vi^ denotes the Christoffel symbol given by 

r;^. = Ig^" {dpga. + d,gp, - d,gp,) (25) 

and r is an affine parameter such that for time-hke geodesies dr^ = g^jjdx'^dx'^ corresponds to proper time. 
The geodesic Lagrangian £g for a point particle in the space-time p^ reads 



1 dx'' dx" _ 11 



-fm'^'M^'ri^hr-H^'' 



(26) 



where (5 = for massless particles and 5 — 1 for massive particles, respectively. 

The constants of motion are the energy E and the angular momentum L of the particle (absolute value and 
direction). To fix the direction, we choose 9 ~ tt/2. Hence the motion will occur in the equatorial plane. The 
conserved quantities then are 

^ „dt I „ dip 

Using these constants of motion we can rewrite ((26l) as follows 



where VcB{r) is the effective potential 



'£'f^l(E'-VMr)). (28) 



VMr) = f{s + {^y (29) 

We then find 

f dr Lf 

The solution to this will give r(ip) for a test particle. Since / and I are only given numerically, r{ip) can only be 
determined numerically. 

IV. NUMERICAL RESULTS 

We have solved the equations of motion ([TB1) - (|18I) numerically using a Newton-Raphson method with an adaptive 
grid scheme [S^]. The relative errors of the numerical integration are on the order of 10~^° — 10"^'^. The geodesic 
equation is integrated numerically in Fortran using the Fortran Subroutines for Mathematical Applications-IMSL 
MATH/LIBRARY. The accuracy of the integration as estimated from application of the method for the integration 
of geodesies in the Schwarzschild space-time is on the order of 10~^. 

A. The metric functions 

In order to be able to compare the space-time of a boson star and that of a Schwarzschild black hole, we have first 
investigated the difference between the metric functions / and I. 



K 


CJ 


Q 


M 


<^(0) 




0.1 
0.1 
0.1 
0.1 
0.1 


0.9260 
0.8820 

0.6800 (1) 
0.6525 

0.6800 (2) 


101.8523 
100.8126 
119.6961 
108.3567 
86.4696 


100.6465 
99.7029 
114.3021 
106.8518 
92.3413 


0.4945 
0.8463 
2.1891 
2.7331 
3.2230 


60.4148 
59.9558 
66.3377 
63.0380 
56.3793 


1.0 
1.0 
1.0 


0.8600 
0.7620 
0.8500 


15.5600 
11.1205 
6.8865 


15.1388 
11.6461 
8.3518 


0.2694 
0.7175 
1.2718 


80.9035 
67.1487 
52.0577 



TABLE I: The values of cj as indicated by the black dots in Fig[2] together with the corresponding values of the charge Q, mass M and the value 
0(0). The (1) and (2) denote the solutions on the 1. and 2. branch, respectively. See also [ 35 || for further details using the potential (fsl. 

For that, we have to remind ourselves of how the Schwarzschild metric looks like in isotropic coordinates which we 
use in this paper. This is given by (J12p with 



fir) 




lir) 






(31) 



where M is the ADM mass of the solution and the isotropic coordinate r is related to the Schwarzschild like coordinate 
f hy f = r{l + M/(2r))^. In FigUwe show the metric functions /(r) and l{r) for a boson star with given mass M 
and in comparison to that the metric functions for a Schwarzschild solution with the same M. 
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FIG. 1: Wc show the metric functions f(r) (left) and l(r) (right) for a boson star with given mass M ^ 10 for k — 1.0. Also given are the metric 
functions /(r) and l(r) for a Schwarzschild solution with the same value of M. Note that the Schwarzschild solution has 

fi'^h) — K^h) — at the event horizon r^ — M/2. 

Clearly, the metric functions agree at large r, while for small r they deviate strongly. In particular the metric 
functions /(r) and l{r) of the Schwarzschild solution tend to zero at the event horizon r = r^ — M/2, while for boson 
stars f{r) and l{r) stay perfectly regular for all r £ [0 : oo[. This fact will play an important role in the behaviour 
of test particles at small r, while our results already indicate here that the behaviour at large r will be similar when 
comparing boson star space-times with the Schwarzschild space-time. 



B. Mass, charge and radius 



It has been observed before that boson stars exist only in a limited range of the frequency uj. At the maximal value 
of UJ both the mass and the charge of the boson stars tend to zero, while at the minimal value of oj a second branch 
of solutions starts to form and finally ends at the critical value of w after a number of spirals. In the following, we 
want to discuss the motion of test particles in boson star space-times. Here we concentrate on the boson stars at 
some particular values of the frequency uj. These are indicated in Fig 121 where we plot the mass M as function of lo 
for K = 0.1 and k = 1.0, respectively. The exact numerical values of the mass M and charge Q are given in Table HI 
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FIG. 2: The value of the mass M is shown as function of uj for two different values of k. The black dots on the curves indicate the values of i 

that we have investigated in more detail and that are given in Table ^ 



In the following, it will also be important to understand what the radius of the boson stars and in particular what 
the extend of the orbits around the boson star are. We would like to discuss this briefly in the following. Since the 
scalar field has an exponential decay in our model, it is difficult to define a definite radius of the boson stars discussed 
here. This is e.g. different in models with a T^-shaped potential [28|, where a definite radius of the object exists. 
Hence, different notions of the boson star radius can be evoked. Either, we can define the radius of the boson star 



(m) 
BS 



"^ (which in our rescaled variables would be equal to unity) or 

(99) 



as the inverse of the scalar boson star mass r 

we could define the radius of the boson star as the radius rjj'g "* inside which 99% of the mass of the star is contained. 
It remains to be shown that stars consisting out of scalar fields can form. Recent results suggest that a fundamental 
scalar boson with mass on the order of to « 125 GeV could exist in nature [3J]. If boson stars would consist of this 

scalar boson, their typical radii would be (in dimensionful quantities [33]) f\^ ~ 10^^^ m. If the star would consist 

of pions with masses of roughly 140 MeV [l^, the typical radius would be fg™ « 10^^^ m. One could also think 
of much lighter scalar particles, e.g. the axion with roughly 10~^ eV and the dilaton with 10""'^° eV which would 
correspond to much larger radii fg™ « 10^^ m and fg™ « 10'^ m, respectively. When computing the radius r^g we 
find that this is roughly two orders of magnitude larger (see Table U for some values). On the other hand, the masses 
of these boson stars can be quite high. The mass M is measured in units of M'^^/m. We find that the maximal value 

of the mass Aimax is well approximated by 



M„ 



7 



K 



47r77|usY"^ ' 



(32) 



which should be compared to the typical mass of a neutron star given by M^ylmf^^^^^^^ with rnncutron ~ IGeV. 



For m sa 10^ GeV this would correspond to masses of lO^^GeV ^ lO^^kg - and this inside a radius of r 



(m) 
BS 



lo- 



rn (or two orders of magnitude higher if we use r^g ). This is an extremely dense object, much denser than a typical 
neutron star. However, these boson stars are stable towards gravitational collapse and due to their extremely high 
density have been considered as alternatives to supermassive black holes [29[. The results for radii and masses in 
dependence on the scalar particle mass are summarized in Table |TT] for a dimensionless value of the boson star mass 
of Af = O(102). 

The relation between the mass M and the charge Q as given in Table U also allows us to make conclusions about 
the stability of the boson stars. Since Q corresponds to the particle number, i.e. the number of massive scalar bosons 
and these have mass m = 1 (in our rescaled variables), an assembly of Q massive scalar particles would have mass 
M = mQ = Q. Hence, as long as the mass and charge of the boson star fulfill M < Q they are stable with respect to 
the decay into Q individual scalar bosons. As such the boson stars on the first branch of solutions are always stable. 



Particle 


m 


M 




Higgs 


125 GeV 


10" kg 


10-1* m 


Pion 


140 MeV 


10" kg 


10-15 m 


Axion 


10-5 eV 


10^^ kg 


10-2 m 


Dilaton 


10-1° eV 


10^2 kg 


10^ m 



TABLE II: The approximate values of the mass M and r^g of the boson star in dimensionful units for different scalar particle with mass m and 
dimensionlcss boson star mass Af — 0(10 ). See also I33II for further details using the potential i fsl . 



while they become unstable on parts of the second branch. 



C. The effective potential and possible orbits 



In all cases, we need to require that E^ > Vcs in order to have solutions to (150]) . At the intersection points E'^ — 14ff 
we have ^ = 0. These are the turning points of the motion. Hence, the effective potential ((29|) can already give 
information about the type of orbits possible. In the following, we want to discuss which type of orbits are possible in 
the space-time of a non-compact boson star. We distinguish different type of orbits as follows: (a) a bound orbit is an 
orbit that possesses a minimal and a maximal finite radius, (b) an escape orbit is an orbit that possesses a minimal 
finite radius, but extends to infinity, (c) a terminating escape orbit is an orbit whose minimal radius is equal to zero 
and on which the particle comes from infinity to then end at the center of the boson star and (d) a terminating orbit 
is an orbit with finite maximal radius and minimal radius equal to zero. The possible orbits for massive test particles 
are summarized in Table IIIIl 

Note that in the following we will concentrate on orbits appearing at reasonably small r. Since the space-time of 
the boson star tends asymptotically to the Schwarzschild space-time, we have in addition all possibilities of orbits 
available there at large r. We will not discuss these in detail, but refer the reader to the literature (see e.g. [36ll37|). 

In Figl3]we show the effective potential VeS for massive and massless test particles, respectively, for a typical boson 
star space-time with k = 0.1 and tu ss 0.68 for the 1. branch as well as for the 2. branch of solutions. 

Again we would like to compare this to the situation in a Schwarzschild space-time. First note that f/l = (1 + 
Af/2r)-^ in the Schwarzschild case such that for r ^ r^ — M/2 the effective potential vanishes since / vanishes. 
As soon as the particle has crossed this horizon we cannot observe it anymore and it will eventually end up in the 
physical singularity at r = 0. For the boson star space-time on the other hand / and I are well behaved and in 
particular never become zero in r S [0 : oo[. As such the angular momentum term I? jr^ dominates the effective 
potential at r — )■ 0. A particle with L ^ can never reach r = no matter how large its energy is. Now depending 
on the massive test particle's energy E and angular momentum L 7^ it will either move on a bound orbit or on 
an escape orbit. Clearly, we can have two intersection points between E^ and V^off for massive particles as soon as 
the effective potential possesses a local minimum with value smaller than the asymptotic value YcsiT = 00) = 1 (see 
Fig 3(a) and Fig 3(c) I. These two turning points at finite values of r correspond to the nfinimal radius Tmin and the 
maximal radius Tmax of the bound orbit. For both values of k that we have studied we observe that the further we 
move on the branches starting at a; = 1, the lower is the value of the minimum of the effective potential. This can 
be clearly seen when comparing the effective potential for a massive test particle with fixed L^ for uj — 0.68 on the 1. 
branch with that on the 2. branch. This means that we can have bound orbits for smaller values of E for a given L 
when moving along the branches. For massive test particles with L ^ 0, but E^ large only one intersection point is 
possible. These particles will then move on escape orbits. 

On the other hand a radially moving massive particle with L — can reach r = if £'^ > V^ff (r ~ 0) which gives 
the condition E'^ — /(O) > 0. Hence, massive test particles without angular momentum and E^ < Vcs{r — 00) = 1 
move on terminating orbits, while those with E^ > 1 move on terminating escape orbits. See also Table IIIIl for the 
possible orbits of massive test particles. 

In Fig |3(b)| and Fig |3(d)| we show the effective potentials for a massless test particle for different values of L in the 
space-time of a boson star with to — 0.68. We observe that the potential is monotonically decreasing and that no 
extrema of the potential exist for L^ > 0, while 14ff(r) = for massless test particles without angular momentum, 
i.e. L^ = 0. Hence, we can only have escape orbits for L^ > and terminating escape orbits for L^ = 0. This is 
summarized in Table ITVl The difference between the 1. branch solutions and 2. branch solutions does not seem to be 
very big for massless test particles. 

For r ^ 1 the metric functions of both the boson star space-time and of the Schwarzschild case approach each 
other. Hence, the effective potential which depends only on the metric functions looks the same in both cases and in 
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FIG. 3: We show the eifcctive potential l |29| l for massive test particles (left) and massless test particles (right), respectively, for different values 
of L". Here k — 0.1 and cj ps 0.68. The plots on the left are for the boson star solutions of the 1. branch, the ones on the right for the 

boson star solutions of the 2. branch. 



No. of turning points 


L^ 


range of r 


orbit 













L" = U 




tcrmmatmg escape orbit 


1 










L" = U 




terminating orbit 


2 










L" > U 




bound orbit 


1 










L" > U 




escape orbit 



TABLE III; Types of orbits possible for massive test particles moving in the space-time of a boson star. 



particular Vcff{r —>■ oo) — > 1 for both the boson star and the Schwarzschild space-time. 



D. Examples of orbits 



In the following, we want to show examples of the possible orbits mentioned above. Bound orbits of massive 
test particles are of particular interest for two reasons: (a) the perihelion shift of the orbits can be compared to 
observational data and (b) black holes are believed to form accretion discs around their centers and since boson stars 
have been considered as an alternative to supermassive black holes something similar should also be possible in that 
case. On the other hand, escape orbits of massless test particles, i.e. in particular photons can be used to compare 
with observational data of the light deflection. 
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No. of turning points L^ 


range of r 


orbit 


1 


U L" =ii 




terminating escape orbit | 


- 1 


i L- >0 




escape orDit | 



TABLE IV: Types of orbits possible for massless test particles moving in the space-time of a boson star. 



1. Massive test particles 

In Fig|4]we show a bound orbit and an escape orbit of a massive test particle with L^ = 0.5 in the space-time of a 
non-compact boson star with k — 0.1 and the values oi uj = 0.68. We compare the motion of massive test particles in 
the space-time of the boson star solutions on the 1. and 2. branch, respectively. 



E'=0.95, L'=0.5 



E'=0.95, L'=0.5 





-10 -5 5 10 



(a) 



bound orbit (1. branch) 



(b) 



bound orbit (2. branch) 



E =1.1, L =0.5 



10 20 30 



E =1.1, L =0.5 



10 20 30 



(c) 



escape orbit (1. branch) 



(d) 



escape orbit (2. branch) 



FIG. 4: We show bound as well as escape orbits for massive test particles with angular momentum L — 0.5 in the space-time of a boson star 
with K — 0.1 and uj — 0.68. For bound orbits the particles have E — 0.95, while for escape orbits we have chosen E — 1.1. The figures 
on the left and right, respectively show the orbits in the space-time of the boson star solution of the 1. branch and 2. branch. In all cases 
the center of the boson star is located at the origin of the coordinate system. 



For the bound orbits we have chosen E^ = 0.95 < V^fi(r = oo), while for the escape orbits we have chosen 
E'^ = 1.1 > Ves{r = oo). Though the frequency is the same in both cases, we see that there are differences in the 
particle motion. The bound orbits differ clearly in their perihelion shift. The perihelion shift in the space-time of 
the boson star of the 2. branch is much larger than for that of the boson star on the 1. branch (see also discussion 
below). The maximal radius of the orbit in the case of the 1. branch on the other hand is slightly bigger (see Table 
|V]for exact numerical values). Moreover, we observe that the further we move on the branches, the lower is the value 
of Tinin, while for k = 1.0 the value of rmax decreases when moving along the branches. For k = 0.1 there seems to be 
a connection between the value of rmax and the value of M with rmax being largest for the maximal possible value of 
M (see Fig 12]). 

The escape orbits also differ in the sense that in the case of the 1. branch the test particles approach the boson star 
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E^ 


L^ 


K. 


U) 


7'miii 


T'max 


e 


0.97 
0.95 
0.95 
0.95 
0.95 


1.0 
1.5 
0.5 
1.5 
0.5 


0.1 
0.1 
0.1 
0.1 
0.1 


0.9260 
0.8820 

0.6800 (1) 
0.6525 

0.6800 (2) 


1.7688 
1.8418 
0.2376 
0.3106 
0.0823 


19.5923 
11.3479 
14.5724 
12.3623 
11.8867 


0.9959 
0.9867 
0.9998 
0.9997 
0.9999 


0.97 
0.97 
0.95 


1.0 
1.0 
1.0 


1.0 
1.0 
1.0 


0.8600 
0.7620 
0.8500 


0.8302 
0.2520 
0.5855 


27.6181 
22.3061 
9.8599 


0.9995 
0.9999 
0.9982 



TABLE V: The values of r 



"max as well as of the eccentricity e :— A/C^J^ax ~ ^miii)/^max °^ bound orbits with given E and L in boson star 

space-times with k and cj. The L value is chosen close to the upper boundary up to which bound orbits are possible such that the 

eccentricity of the orbits is close to its maximal possible value. 



and get deflected into a direction nearly perpendicular to the original one. In the case of the 2. branch the particle 
encircles the boson star and gets deflected into a direction similar to the one that it came from. 

a. Innermost stable circular orbits (ISCOs) If the effective potential possesses a local minimum and we choose 
E"^ such that it is equal to the value of the effective potential at this minimum the orbits are circular orbits, so-called 
"innermost stable circular orbits" (ISCOs). Examples of these orbits are shown in FiglS]for different values of w and 
L^ . Note that as soon as we have fixed the L^ the ISCO exists only for one specific value of E^ . 

We observe that the further we move on the branches, the smaller is the radius of the ISCO of a particle of given 
angular momentum. This can be seen for a test particle with L^ = 1.5 when comparing Fig 5(b) and Fig 5(e) as well 
as for a particle with L^ = 0.5 when comparing Fig 5(c) and Fig 5(d) Apparently, the extend of the orbit is not only 
connected to the mass of the boson star since the ISCO of a L'^ — 1.5 particle has larger radius in the space-time 
of a boson star with oj — 0.8820, M = 99.7029 as compared to that in a boson star space-time with lo = 0.6525, 
M — 106.8518, while the ISCO for a L^ = 0.5 particle has a larger radius in a boson star space-time with cu — 0.6800, 
M — 114.3021 as compared to that in a boson star space-time with cj = 0.6800, M = 92.3413. The message that one 
should take away here is that the non-trivial scalar field seems to play an important role in the particle motion and 
that no direct relation between the mass of the object and the ISCO exists as one would expect from a Newtonian 
treatment of the orbits. 



2. Massless test particles 

In Figiniwe show the escape orbits of massless test particles with E'^ =0.1, L^ = 1.0 in the space-time of the boson 
stars given in Table HI Depending on the choice of cj and k wc find that the particles get deflected into a direction 
between roughly 7r/4 and 7r/2 and that for k large enough the minimal radius of the escape orbit lies inside the 
Schwarzschild radius of the corresponding black hole with the same mass. We do not plot terminating escape orbits 
of massless test particles here since they would just correspond to a line from infinity to r = 0. 



E. Observations of boson stars 



Next to the possibility that boson stars could act as a toy model for neutron stars, boson stars have been suggested 
as alternatives to supermassive black holes residing in the center of galaxies [29|. As mentioned above, objects with 
a well-defined surface as alternatives to the supermassive black hole at the center of our own Milky Way have been 
ruled out [3l[. We have, however argued above that our boson stars are non-compact in the sense that there is no 



well-defined surface beyond which the energy density and pressure, respectively, of the boson star strictly vanishes. 



1. Boson star at the center of our galaxy 



The center of our Milky Way possesses a very compact astronomical radio source named Sagittarius A* . It is 
believed to be a supermassive black hole as observations of the motion of stars around the galactic center suggest 
that the mass of the central object is on the order of M w IQ^'^kg which corresponds to roughly 4 • 10® solar masses 
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E =0.884, L =1.0 



E =0.891, L =1.5 





-2 2 



(a) uj ^ 0.9260, L^ ^ 1, E^ ^ 0.884 (b) uj ^ 0.8820, L^ ^ 1.5, B^ ^ 0.891 



E''=0.489, L''=0.5 



E''=0.405, L''=0.5 





(c) cj ^ 0.6800 (1. branch), L^ ^ 0.5, £;^ ^ 0.489 (d) cj ^ 0.6800 (2. branch), L^ ^ 0.5, E^ ^ 0.404 



E''=0.742, 1^^=1.5 




(e) UJ ^ 0.6525, L^ ^ 1.5, £;^ ^ 0.742 



FIG. 5: Wc show nearly circular, so-called "innermost stable circular orbits" (ISCOs) for massive test particles with given angular momentum 

L in the space-time of a boson star with k — 0.1 and different values of to. 
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E'^=0.1,L'^=1.0 



E''=0.1,L'^=1.0 



-20 20 
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(a) 



u) = 0.9260, K = 0.1 



(b) 



u) = 0.8820, K = 0.1 



E-^-O.I , L-^.l .0 



£■^-0.1 , L-^.l .0 



-20 20 




10 20 30 



(c) 



bJ = 0.6800 (1), K = 0.1 



(d) 



ui = 0.6525, K = 0.1 



-20 20 



(e) 



(.J = 0.6800 (2), K = 0.1 



(f) 



ui = 0.8600, K = 1.0 



E^-O.l , L'.l .0 




E^.0.1 , L^.1 .0 



-20 -10 10 20 



(g) 



ui = 0.76200, K = 1.0 



(h) 



uj = 0.8500, K = 1.0 



FIG. 6: Wc show escape orbits of niasslcss test particle with E — 0.1 and L — 1.0 in the space-time of a boson star with k — 0.1 and k — 1.0, 
respectively for the boson stars given in Tabic ^ Also indicated is the radius of the event horizon of a Schwarzschild black hole with the 
same mass as that of the boson star (red circle) if this radius is larger than the minimal radius of the escape orbit, which is the case only 

for K = 1.0. 
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[3a, 1231, which is contained within a sphere of radius roughly 2 • lO^^m. Comparing with Table HIl we observe that we 
can have boson stars that can be as dense as the observed object. Furthermore, since the orbits of a number of stars 
have been determined, we can compare the observational results with ours. Unfortunately, the perihelion shift of the 
orbits of stars around the galactic center has not yet been determined (further observations are necessary here), but 
the eccentricity e := \/(?'niax ~ '"min)/''max ^^ the orbits has been calculated from the observations. In Table|V]we give 
some values of e for orbits in boson star space-times. Note that these are the nearly maximally available eccentricities, 
but that we can also have orbits which are nearly circular (see discussion above on ISCOs). Hence e > up to e < 1 
is possible. Interestingly, some of the orbits of stars circling around the center of our own Milky Way ha ve q uite large 
eccentricities |40l - |42| . e.g. the orbit of the star called 5*14 could have eccentricity up to e « 0.974 ±0.016 [40] • Clearly, 
we are able to describe these large eccentricities with our model. 

Another interesting observation has been done related to the emission of short wavelength radio waves from Sagit- 
tarius A* . These are emitted due to accretion of surrounding matter onto the supermassive object in the center of 
our Milky Way. In I43| radio waves with wavelength of 1.3 mm have been observed in order to understand the region 
close to the supermassive object in which gravitational effects should be very strong. Observing Sagittarius A* at 
wavelength 1.3 mm, the authors of [43| concluded that the size of this object is smaller than the expected apparent 
event horizon size of the assumed black hole. It was suggested that this might be due to the fact that the biggest 
fraction of emission of radio waves is not due to the black hole at the center, but rather related to the accretion in 
the vicinity of the supermassive object. 

If we would, however, give up the idea of a black hole in the center of our galaxy and use a non-compact boson 
star instead particles could be emitted from radii smaller than the apparent event horizon of the black hole since the 
space-time is globally regular. This is demonstrated in Figl?! where we show bound and escape orbits of massive test 
particles in the space-time of a boson star with k = 1.0 and different values of uj. For comparison wc also show a 
circle with radius equal to the Schwarzschild radius of a black hole with the same mass, which in our coordinates is 
rh = M/2. It is apparent that the minimal radius of these orbits (both bound as well as escape) is smaller than the 
event horizon, while the maximal radius of the bound orbits is larger. If we would hence have stricter bounds on 
the mass (and hence the corresponding Schwarzschild radius) of the supermassive object in the center of our galaxy 
and would have clear signals from particles coming from "inside" the event horizon a possible explanation might be a 
non-compact boson star. Future observations at wavelength even smaller than 1.3 mm should be able to show what 
the space-time of the supermassive object at the center of our galaxy looks like. 



2. Perihelion shift 

The motion of a number of stars orbiting the center of our Milky Way have been determined [40l - |43 |. Since the 
orbital periods To of these stars are quite large (between To = 11.5±0.3 years for SO-102 [42| and To = 94.1±9.0 years 
for SI [4l|) it will need longer observations to determine the perihelion shift of these objects. But this is certainly one 
of the interesting observables. Hence, we discuss the possible perihelion shifts of the orbits of massive test particles 
in the space-time of non-compact boson stars here. These might be later used to compare with observations. 

The perihelion shift of a bound orbit of a massive test particle is given by 

^miii 

with Ves given by ([3]), where we set (5 = 1. 

The values for orbits with nearly maximal eccentricity are given in Table IVll We observe that the perihelion shift 
in all cases is quite large and can even become negative. This latter fact would be something that is not possible in 
a Schwarzschild black hole space-time and would allow us to distinguish black holes from boson stars. 



3. Light deflection 

Since the gravitational field close to a supermassive object is very strong it is important to understand light deflection 
in order to be able to interpret observational data correctly. Due to the strong bending of light rays the apparent 
event horizon of the proposed black hole would be twice as big as the actual event horizon. All of these things have 
to be taken into account when analyzing observational data. 
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(e) 



bound orbit, i.j = 0.8500 



(f) 



escape orbit, lj — 0.8500 



FIG. 7: We show bound and escape orbits of massive test particle in the space-time of a boson star with k — 1.0 for different values of lj. Also 
indicated is the Schwarzschild radius r/^ — AI/2 (red circle), i.e. the radius of the event horizon of a black hole with the mass M equal to 

that of the respective boson star (see also Table QJ- 



UJ 


K 


E^ 


L^ 


(Av9),=i 


0.9260 


0.1 


0.97 


1.0 


-0.51227r 


0.8820 


1.0 


0.95 


1.5 


-0. 144871 


0.6800 (1) 


0.1 


0.95 


0.5 


-0.22847r 


0.6525 


0.1 


0.95 


1.5 


1.35667r 


0.6800 (2) 


0.1 


0.95 


0.5 


0.49337r 


0.8600 


1.0 


0.97 


1.0 


-0.416l7r 


0.7620 


1.0 


0.97 


1.0 


0.48527r 


0.8500 


1.0 


0.95 


1.0 


1.23867r 



TABLE VI: The values of the perihelion shift (Aip)s=i for the boson star space-times given in Table[T]and for the orbits with nearly maximal 

eccentricity (see also Table \V\ . 
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OJ 


K 


E^ 


L^ 


(A(^)i=o 


0.9260 


0.1 


0.1 


1.0 


0.16197r 


0.8820 


0.1 


0.1 


1.0 


0.19367r 


0.6800 (1) 


0.1 


0.1 


1.0 


0.28297r 


0.6525 


0.1 


0.1 


1.0 


0.62287r 


0.6800 (2) 


0.1 


0.1 


1.0 


0. 167871 


0.8600 


1.0 


0.1 


1.0 


0.35687r 


0.7620 


1.0 


0.1 


1.0 


0. 303971 


0.8500 


1.0 


0.1 


1.0 


0.16367r 



TABLE VII: The values of the light deflection (Aip)5=o for the boson star space-times given in Tableland E — 0.1 and L — 1.0. 

The light deflection in the space-time of a boson star is given by 

oo 

dr Lf 



(Av3)5=o = 2 



^E^ - V^sir) r^Vl 



(34) 



with T4ff given by ([3]), where we set (5 = 0. 

In Table rvm we give the light deflection in the boson star space-times indicated in Table H] for massless test particles 
with E^ — 0.1 and L^ — 1.0. We observe that the light deflection increases when moving along the branch for 
K — 0.1 up to roughly the minimal possible value of w « 0.6525. On the second branch, the light deflection decreases 
again. Similar to what was stated above for massive test particles, there seems to be no direct correlation between 
the mass of the boson star and the value of light deflection and hence the presence of the scalar field seems to have 
a non-trivial effect. For k = 1.0 we find that the light deflection decreases when moving along the branches. This 
again indicates that a semi-analytical treatment of the space-time of boson stars does not lead to the proper results 
and that numerical techniques are vital to make predictions about observables in such kind of space-times. 



CONCLUSIONS AND OUTLOOK 



In this paper, we have studied the geodesic motion of test particles in boson star space-time and have discussed 
possibilities to differ them from the geodesies around a Schwarzschild black hole. Our work has many interesting 
applications. As a first example let us mention the accretion of matter onto massive compact objects such as neutron 
stars. Neutron stars are very dense objects that are typically modeled using an appropriate equation of state. Boson 
stars could act as a toy model for these objects. In the following, one could use our techniques to model accretion 
flow onto a boson star in a simple ballistic a ppr oach. This was done before in black hole space-times where analytic 
solutions to the geodesic equation exist [4J, |4^ . This technique could be applied as well using the numerical data 
of a boson star space-time. The accretion flow would also be interesting from the point of view of interpreting the 
boson star as an alternative to a supermassive black hole residing in the center of our Milky Way. The accretion of 
matter leads to the emission of radio waves and these can be used to make predictions about the properties of the 
supermassive object at the center of our galaxy. Gravitational wave signals would be another possibility. These have 
already been discussed in |30j . 

It would also be interesting to see how test particles move in radially excited boson star space-times as well as in the 
space-time of a rotating boson star. For once this is interesting since it is believed that the flat space-time counterparts 
of boson stars, so-called Q-balls could originally form in an excited state Q. Furthermore, due to the interaction with 
surrounding matter, e.g. accretion, boson stars might be excited. If boson stars should act as toy model for neutron 
stars it would be surely also of interest to study the rotating counterparts since basically all observed neutron stars 
seem to rotate, in fact quite quickly as the example of a radio pulsar (assumed to be a neutron star) spinning at 716 
Hz shows [46|. It would then be interesting to study the motion of test particles in these space-times. 

Finally, in order to model accretion flow onto dense astrophysical objects more realistically one should consider 
charged test particles (that would make up the plasma). Hence the equation of motion for the particle should be 
modified to include the charge of the particle. Studies of charged test particles in some black hole space-times have 
been done previously [43| . In order to have interactions of the test particles with the boson star the latter should be 
charged as well. This can be done by promoting the global U(l) symmetry of our model to a local U(l) and has been 
considered for a number of boson stars H&\ . The corresponding equation of motion for the test particle would then 
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have to be solved numerically. 
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